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Abstract 

Background: Information processing in the brain requires large amounts 
of metabolic energy, the spatial distribution of which is highly hetero- 
geneous reflecting complex activity patterns in the mammalian brain. 
Results: Here, it is found based on empirical data that, despite this 
heterogeneity, the volume-specific cerebral glucose metabolic rate of 
many different brain structures scales with brain volume with almost 
the same exponent around —0.15. The exception is white matter, the 
metabolism of which seems to scale with a standard specific exponent 
— 1/4. The scaling exponents for the total oxygen and glucose con- 
sumptions in the brain in relation to its volume are identical and equal 
to 0.86±0.03, which is significantly larger than the exponents 3/4 and 
2/3 suggested for whole body basal metabolism on body mass. 
Conclusions: These findings show explicitly that in mammals (i) 
volume-specific scaling exponents of the cerebral energy expenditure 
in different brain parts are approximately constant (except brain stem 
structures), and (ii) the total cerebral metabolic exponent against 
brain volume is greater than the much-cited Kleiber's 3/4 exponent. 

The neurophysiological factors that might account for the regional 
uniformity of the exponents and for the excessive scaling of the total 
brain metabolism are discussed, a^ong with the relationship between 
brain metabolic scaling and computation. 
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Background 

The brain is one of the most expensive tissues in the body [1, 2, 3], as it uses 
large amounts of metabohc energy for information processing [4, 5, 6, 7]. Because of 
this, neural codes are constrained not only by a combination of structural and func- 
tional requirements [8, 9, 10, 11, 12, 13, 14, 15, 16, 17], but also by energetic demands 
[6, 7, 18, 19, 20, 21]. In general, it is has been observed that an elevated synaptic sig- 
naling between neurons leads to more energy consumed [22], which is used in imaging 
experiments [23] . Although some theoretical progress has been made in quantifying con- 
tributions of different neurophysiological processes to the total metabolic expenditure 
of a single neuron [20] , the metabolism of large-scale neural circuits has not been inves- 
tigated quantitatively Here, a first step in this direction is made by studying global and 
regional in vivo brain metabolic scaling. There is a long tradition in applying allometric 
scaling to problems in biology [1], in particular to whole body metabolism [1, 24], but 
surprisingly not to cerebral metabolism. The goal of this paper is to find, by collecting 
and analyzing data, scaling metabolic exponents of different parts of the brain, as well 
as its global exponent. It is found that the volume-specific exponents across cerebral 
regions on brain volume are almost identical, approximately —0.15. Consequently, the 
energy consumption of the entire brain tissue scales with brain volume with the expo- 
nent ~ 0.86. The main neurophysiological factors that might cause the increase of the 
latter exponent well above the putative Kleiber's 3/4 scaling exponent characterizing 
whole-body energy expenditure on body mass [24], are identified. The consequences of 
brain metabolic scaling on its information processing capacity in the context of brain 
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design are also discussed. 
Results 

Oxygen and glucose are the main components involved in the production of ATP, 
which is used in cellular energetics [3, 4, 6], and therefore their rates of utilization provide 
useful measures of brain metabolism. There are several mammalian species, spanning 
more than 3 orders of magnitude in brain size, for which in vivo brain metabolic data 
are available (see Materials and Methods). The allometric laws characterizing global 
cerebral metabolism of oxygen and glucose are similar and yield an identical scaling 
exponent 0.86 ± 0.03 against brain volume (Fig. 1). It is important to note that this 
value is significantly larger (p < 0.05) than the exponents 3/4 [1, 24] and 2/3 [25] found 
for whole body mammalian metabolism in relation to body mass. 

The cerebral cortex is a critical part of the brain responsible for integrating sensory 
information, and commanding behavioral and cognitive tasks. Regions of the cerebral 
cortex differ both in molecular detail and in biological function, which is manifested 
in a non-uniform distribution of neuronal activity and energy utilization throughout 
the cortex [5] (and Supplementary Information). However, despite this heterogeneity 
values of the scaling exponents of the regional volume-specific glucose utilization rates 
on brain volume {CMRgic, glucose cerebral metabolic rate per brain region volume) are 
surprisingly homogeneous; they are either exactly equal to —0.15 or are close to this 
value (Fig. 2 and Table 1). Consequently, also the average specific glucose utilization 
rate of the whole cerebral cortex scales with the exponent —0.15 ± 0.03 with brain 
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volume (Fig. 2E), which is equivalent to the exponent 0.85 ± 0.03 for the metabohsm 
of the entire cortical volume; the value close to that for the whole brain (Fig. 1). 

Some subcortical structures of gray matter utilize half of the energy used in the cor- 
tex (e.g., limbic structures in cat and monkey; see Clarke and Sokoloff [5] and Supple- 
mentary Information), and yet almost all of them exhibit a similar scaling homogeneity, 
with metabolic specific exponents also around —0.15 (Fig. 3, and Table 1). Volume- 
specific metabolisms of two brain stem structures - superior coUiculus (involved in visual 
coordination) and inferior coUiculus (involved in auditory processing) - seem to be ex- 
ceptions, since they scale with brain size with the exponent ^ (Table 1). This might 
be caused by their highly variable activities (see Supplementary Information), and we 
do not know how other brain stem structures behave. The high degree of allometric 
uniformity for the most of the subcortical system is even more striking than that for 
the cortex, since subcortical regions are much more diverse in function and biophysical 
properties than cortical areas. For example, thalamus and hippocampus play extremely 
different roles in the brain, the former mediating sensory input to the cortex, while the 
latter implicated in memory processes, and still their scaling exponents and correspond- 
ing confidence intervals are almost identical (Figs. 3A,B, and Table 1). 

Metabolism of white matter is about three-fold lower than that of gray matter 
[4, 5], and the results in Table 1 indicate that also their scahng exponents might differ. 
Specific glucose metabohsms of the two main structures of white matter, corpus callo- 
sum and internal capsule, scale similarly with brain volume with the average exponent 
around —1/4 (Fig. 4, Table 1). This value resembles the whole body specific metabohc 
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exponent against body mass. 



Discussion 

General discussion. 

The total metabolic exponent 0.86 ± 0.03, found for unanesthetized mammalian 
brains, implies that cerebral energy use increases more steeply with brain size than 
does whole body energy use with body size. This might be a reason why brain size 
increases slower than its body size, with the power 3/4 [26]. However, the metabolic 
exponent 0.86 undermines Martin's [26, 27] well-known argument that the 3/4 scahng 
of brain size reflects allometric isometry between brain metabolic needs and its size 
(exponent 1). Instead, the product of the exponents for brain size on body size and 
brain metabolic rate on brain size gives an exponent for brain metabolic needs on body 
size of about 0.65, which is lower than 0.75, expected from the assumption of isometry. 
Thus, depending on the scaling reference the total brain metabolism scales either above 
the Kleiber's 3/4 power if in relation to brain mass, or below the Kleiber's power if in 
relation to body mass. 

The discovered uniformity, i.e. constancy, of the cerebral specific exponents in gray 
matter suggests a common principle underlying basal metabolism of different brain 
structures, which might be associated with the homogeneity of synaptic density through- 
out the gray matter [28, 29]; see below. 

In analyzing comparative allometry some authors use phylogenetic approaches in or- 
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der to include dependencies in data sets [30, 31]. However, these sophisticated methods 
require as a prerequisite the knowledge of a phylogenetic tree and associated branching 
parameters for species of interest. This is not a trivial matter to do, and therefore not 
applied in the present analysis. Also, because correlations in most scaling plots (Figs. 
1-4, Table 1) are generally very high, it is likely that taking phylogeny into account 
would not alter the empirical exponents. 

Key factors in the cerebral metabolic rate. 

Which factors in the gray matter might account for its total metabolic exponent 
on brain volume greater than 3/4? The likely candidate is the density of glial cells, 
which provide metabolic support for neurons (including synapses) [32, 33]. Two recent 
studies [34, 35] show that the number of glia per neuron increases for larger brains, 
suggesting that neurons become more energy expensive with increasing brain size. In 
particular, the total number of glia in the cerebral cortex of rodents scales with brain 
size with the exponent 0.89 (Fig. 2B in [34]), which is close to the empirical metabolic 
exponent 0.86. It is likely that glia number simply follows energetic demands of neurons, 
especially their spatial expansion with increasing brain size. Thus, the details of this 
expansion may provide clue about which neural factor are important. 

It has been estimated that neural metabolism is dominated by Na+/K+-ATPase ion 
transport [6, 36, 37, 38]. The bulk of this comes from active ion fluxes at synapses 
and along axon (propagation of axon potentials) if neural firing rate is sufficiently large 
[20] . The synaptic contribution in a single neuron is proportional to the product of the 
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number of synapses per neuron, firing rate, release probability, and the postsynaptic 
charge. The active axon contribution is proportional to its surface area and firing rate. 
When neurons do not fire action potentials they also consume energy, because of the 
passive Na"'" and ion flow that electro-chemical Na''"/K"'" pump must remove to 
maintain their gradients across the membrane [39]. This resting potential contribution 
is proportional to the total neuron's surface area. To obtain the total cerebral energy 
consumption one has to multiply all these three additive contributions by the total 
number of neurons in the gray matter. 

It seems that from all these three neural contributions the most dramatic for the 
scaling exponent of the total brain metabolism on brain volume is the synaptic contri- 
bution. This contribution is proportional to the total number of synapses in the gray 
matter, which scales with the gray matter volume, and thus brain volume [40, 41], with 
the exponent 1. This follows from regional- and scale- invariance of synaptic density 
[28, 29]. The regional homogeneity of synaptic density correlates with the discovered 
homogeneity of the regional volume-specific cerebral metabolic scaling (Table 1). More- 
over, if remaining factors comprising the synaptic contribution, i.e., firing rate, release 
probability, and postsynaptic charge were brain size independent, then the synaptically 
driven total brain metabolism would scale with brain volume with the exponent 1. Since 
this exponent is, in fact, between 3/4 and 1 (Fig. 1), it implies that either of these 3 
factors (or all of them) decrease with brain size. It could be hypothesized that because 
synaptic sizes and their basic molecular machinery are similar among mammals of dif- 
ferent sizes [28, 39], the postsynaptic component might be roughly the same among 
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different species. Tfiis suggests tfiat tfie product of the firing rate and release probabil- 
ity should decrease as brains increase in size, with a power of about —0.15, which is in 
accord with low firing rates in humans estimated based on their basal cerebral metabolic 
rate [21, 42]. This also suggests that the number of active synapses in the background 
state decreases for bigger brains. That firing rate should decrease with brain (body) 
size is also consistent with allometric data of avian sensory neurons firing rates [43]. 

The remaining two contributions affecting metabolic rate, the active axons and 
maintaining the resting potentials, are both proportional to the product of the number 
of neurons in the gray matter and axonal surface area (assuming that axon surface 
area constitutes the majority of the neuron's area, especially for bigger brains). Addi- 
tionally, the active axon contribution is proportional to the average firing rate. Given 
the above indications that the firing rate likely decreases with brain volume the active 
axon contribution becomes less important for the total metabolic exponent as brains 
increase in size. Thus, we focus only on the resting potential contribution. Since the 
product of the number of neurons and axonal surface area is proportional to the ratio 
of the volume of gray matter to axon diameter (due to the empirical fact that volumes 
of intracortical axons and gray matter are proportional across species [28]), the resting 
potential contribution produces the total metabolic exponent that is determined by the 
scaling exponent of the axon diameter against brain volume. The bigger the axonal 
exponent, the smaller the metabolic exponent. There are some sketchy experimental 
indications that axon diameter indeed increases with brain size [44, 45], e.g., in the 
corpus callosum with the exponent ^ 0.07 [44]. If similar allometry holds for the gray 
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matter axons, the resting potential contribution would yield a total metabolic exponent 
also above 3/4. 

Brain metabolism, computation, and design. 

The facts that the synaptic metabolic contribution decreases with decreasing the 
rate of release probability, and that the active axon and resting potential metabolic 
contributions decrease with increasing axon diameter have interesting functional con- 
sequences. Higher failure rate of synaptic transmission as brains get bigger not only 
saves energy but it also may maximize information transfer to postsynaptic neurons 
[19]. Similarly, increasing axon diameter with brain size accomplishes three functions 
simultaneously: it reduces the specific metabolic rate, increases the number of synapses 
per neuron, and it increases the speed of signal propagation in cortical circuits, which 
is proportional to the square root of the axon diameter [46, 47] . It is difficult to say at 
this stage whether these relationships are accidental or maybe a result of some cerebral 
optimization. 

We can estimate the allometric cost of information processing by finding how the 
amounts of metabolic energy per neuron and per synapse scale with brain size. Since 
the total energy utilized by the entire cerebral cortex (gray matter) scales with its vol- 
ume Vg with the exponent ^ 0.85 (Table 1, Fig. 2E), the cerebral energy per neuron 
is cx Vg'^^/{pnVg) (X \/0-05-o.i7^ increases with brain size, where we used the fact 

that the neural density p„ scales with brain volume with the exponent between —0.32 
and -0.20 (based on data of Haug (1987) [48]; see Supplementary Fig. S2). This in- 
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crease is in accord with the trend of increasing the number of ghal cells per neuron 
in gray matter. The opposite is true for synapses, since the energy per synapse is 
oc Vg-^^ / {psVg) oc V~^-^^, where Ps is the synaptic density (independent of brain size). 
The increase in energy expenditure per neuron with increasing brain size is in con- 
trast to findings in liver cells [49], whose metabolism decreases with body mass. This 
difference reflects the increase of neural size (its wire) and corresponding decrease in 
density with increasing brain volume (sizes of liver cells are roughly constant [50]). The 
decaying trend for synapses implies again that their active fraction decreases as brains 
get bigger. These results have implications for coding and cortical organization. The 
fact that expanding neurons are energetically costly was probably a driving evolution- 
ary force in decreasing their density in larger mammals (Supplementary Fig. S2) and 
adopting sparse neural representations [18, 20, 21]. The latter factor is consistent with 
the idea of functional modularity of the cerebral cortex [40], i.e., that primary infor- 
mation processing takes place in local modules/areas. The sizes of such modules/areas 
seem to follow scaling rules [12, 15, 40], and have been shown to have almost brain size 
independent connectedness [12] as opposed to neural connectedness that decays with 
brain size [11, 12]. 

Metabolism of gray versus white matter. 

The data in Tabic 1 seem to indicate that the white matter and gray matter 
metabolic allometries are different. This finding implies that as brains increase in size, 
the white matter metabolism is less costly than that of gray matter. This is presum- 
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ably beneficial for the total cerebral energetic expenditure, since white matter increases 
disproportionately faster than gray matter [9, 40, 51]. The difference in white and 
gray matter metabolisms may be caused by their apparent neuroanatomical differences, 
since most of the white matter axonal membrane is covered with myelin sheath, which 
prevents ions flow and reduces metabolic cost. 

Relation to metabolism of other tissues. 

Most of the tissues in the body have much lower specific metabolic rate than the 
brain, with the exception of four highly active organs: kidney, liver, heart [1] and 
gut [2]. There exists no reliable data on in vivo allometric metabolic scaling in these 
tissues across different species (see however [52], where allometric exponents are given 
based on only 2-3 species). Allometric in vitro studies of Na^/K^-ATPase in kidney 
[53] and in brain [54] suggest that these two organs might have comparable specific 
metabolic exponents. We can indirectly estimate and compare metabolic exponents for 
active tissues using allometric data on mitochondria size [55] , since its total membrane 
surface area correlates with a baseline oxygen consumption [55] . (Interestingly, the total 
mitochondrial volume in locomotory muscles is proportional to the maximal metabolic 
rate in mammals [56] .) We find that the total mitochondrial surface area in brain scales 
with brain mass with the exponent 0.86, i.e., exactly the same as that in Fig. 1. The 
corresponding scaling exponents for kidney, liver, and heart, against their respective 
organ masses are smaller and closer to the 3/4 exponent: 0.71, 0.74, and 0.81. These 
exponents do not seem to relate directly to the exponents of sizes of these tissues on body 
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size, since masses of kidney and liver increase slower than body mass, while heart mass 
scales isometrically [1]. Thus, it appears that, in general, higher metabolic exponents 
do not necessarily lead to lower mass exponents. 

If these exponents reflect a real difference of metabolic allometry between cerebral 
and non-cerebral tissues, it might be caused by differences in membrane chemical com- 
position and ion pump activities, which is known as the "membrane pacemaker theory of 
metabolism" [38]. For example, it has been shown that heart, hver, kidney, and skele- 
tal muscles display allometric variation in lipid composition but brain does not [57]. 
Other potential factors (some of which might be related to the membrane pacemaker 
hypothesis) affecting differences in the allometries of brain and other tissues include: (i) 
distinct ways energy is utilized in the brain and in other tissues (e.g.. Na+/K'''-ATPase 
dominates energy consumption only in the brain and kidney [3]); (ii) difference in a 
mode of activity (cells outside nervous system exhibit graded electrical activity without 
firing Na+ action potentials); (iii) structural differences between brain and other highly 
active tissues (the size of non-cerebral cells is virtually independent of the body mass 
and the cells lack elongated processes with synapses, e.g., [50]); (iv) the existence of the 
blood-brain barrier that restricts a direct transport of molecules between bloodstream 
and nerve cells [39], which might affect substrate utilization rate and/or neural activity 
[58, 59]. 

Supply-limited models of metabolic scaling. 

In many studies of whole body metabolism the scaling exponent 3/4 was found [1, 
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24] , and it was argued that this value follows from a general model of hierarchical fractal- 
like transport networks [60], or from constrained geometric networks with balanced 
supply and demand [61]. Both of these models are based on the assumption that 
metabolic rates are determined solely by resources supply rates and are independent 
of the cellular energy expenditure. This single-cause assumption has been challenged 
recently [62, 63]. The main arguments against the above supply-limited models are that 
(i) blood flow rate adjusts itself to tissue physiological demands and in resting animals 
is well below its maximal limit, and (ii) cellular metabolic rates decline with increasing 
body size [49] (at least for non-cerebral cells). Because of these facts a "multiple-causes" 
scenario of metabolic scaling has been proposed [62, 63], which in essence argues that 
supply rate is only one of the factors and should be considered together with other 
factors characterizing utilization of cellular energy. The approach taken in this paper 
is similar in spirit, i.e., given that the total cerebral metabolic exponent is 0.86, simple 
supply-limited models are rejected as a possible explanation for brain metabolic scaling. 
Instead, the most energetically expensive cellular factors that are most likely to affect 
the metabolic exponent were identified. In this sense, this approach can also be viewed 
as a multiple-cause model of the cerebral metabolic scaling. 

Conclusions 

Figures 1-3 and Table 1 provide an empirical evidence that the scaling exponents 
describing global and regional brain metabolism are significantly different (p < 0.05) 
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from the much-cited 3/4 exponent, which calls into question the direct applicability of 
supply-limited models [60, 61] to brain metabolism. The exceptions are white matter 
structures, which seem to exhibit "regular" metabohc exponents (Fig. 4). 

The empirical results presented in this paper show striking uniformity of the al- 
lometric metabolic exponents over almost entire gray matter of mammalian brains at 
normal resting conditions, despite anatomical and functional heterogeneity of different 
regions and their different levels of activation. This regional scaling uniformity is sur- 
prising, as activity level could potentially affect the scaling exponent. For example, the 
total metabolic rate of maximally exercised body scales with body mass well above 3/4, 
with an exponent of 0.88 [56, 64]. 

Materials and Methods 

In vivo data of the cerebral oxygen (CMR02) and glucose {CMRgic) utilization 
rates of unanesthetized adult animals at resting conditions were collected from dif- 
ferent sources [65] (see Supplementary Information for details). In those studies the 
measurements of glucose utilization in all species were performed by essentially the 
same method or its modification (in human and baboon), and thus all glucose data 
are directly comparable. There is a small method variability for oxygen data, since the 
same technique was applied to five out of seven species (except cat and dog). However, 
this variability does not affect the scaling exponent (Fig. lA) - it is exactly the same 
even if only single- method mammals are included in the plot. Glucose utilization data 
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represent both global and regional cerebral metabolism. The investigated mammals 
include: Swiss mouse (only glucose data), Sprague-Dawley rat, squirrel (only glucose 
data), rabbit (only glucose data), cat, dog (only oxygen data), macaque monkey, ba- 
boon, sheep, goat (only glucose data), and human. The investigated brain structures 
include: cerebral cortex (visual, prefrontal, frontal, sensorimotor, parietal, temporal, 
cingulate, occipital), thalamus (including lateral geniculate nucleus and medial genicu- 
late nucleus), hypothalamus (and separately mammillary body), cerebellum (including 
cerebellar cortex and dentate nucleus), basal ganglia (caudate, substantia nigra, globus 
pallidus), limbic system (hippocampus, amygdala, septum), brain stem (superior col- 
liculus, inferior coUiculus), and white matter (corpus callosum, internal capsule). In the 
cases when there are more than one data point for a given animal or a brain structure, 
an arithmetic mean of all values was taken. 

AUometric metabohsm of the entire cerebral cortex, presented in Fig. 2E, was 
obtained by computing an arithmetic mean of glucose utilization in 8 cortical areas 
(listed above) for each animal. Glucose utilization of a given area was itself an average 
of values taken from different sources. If data for all 8 areas were not available, averaging 
was performed over lesser number of areas. For consistency, also an alternative method 
of averaging was used: first averaging was performed for a given source data, and second 
among different sources representing the same animal. In this method, because various 
sources differ in the number of cortical areas studied, averaging in many cases was 
performed over significantly different number of areas. However, both methods give 
statistically identical scaling exponents for the cerebral cortex metabolism (see Fig. SI 
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in Supplementary Information). 

AUometric metabolism of the entire brain was obtained by using either direct data 
quoted by authors, or, if not available, by computing an arithmetic mean of glucose 
consumption in all brain structures provided by the authors. For all scaling plots brain 
volumes were taken from Hofman (1988) [41] and Stephan et al (1981) [66], or from the 
source. 
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Specific scaling exponents of the regional cerebral glucose utilization rate CMRgic 



against brain volume. 
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t-test shows that population mean of the exponents in the second column is significantly 



greater than the exponent —0.25 (p = 9-10 ^, df=21 if white matter included; p = 
2 ■ 10~^, df=19 if white matter excluded). 
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Figure Legends 
Fig. 1 

Scaling of the total basal cerebral metabolism with brain volume. The least-square fit 
line for the log-log plot yields: (A) For the total oxygen consumption rate the scaling 
exponent is 0.86 ± 0.04 (y = OMx - 1.02, = 0.989, p < 10"^, = 7), and its 95% 
confidence interval is (0.75, 0.96). (B) For the total glucose utihzation rate an identical 
exponent 0.86 ± 0.03 is found {y = OMx - 0.09, = 0.994, p < 10-^ iV = 10), and 
its 95% confidence interval is (0.80,0.91). 

Fig. 2 

Scaling of the cerebral cortex specific glucose utilization rate, CMRgic, with brain vol- 
ume. The specific metabolic scaling exponent, corresponding to the slope of the regres- 
sion line, has the following values: (A) —0.12 for visual cortex {y = — 0.12x + 0.02); 
(B) —0.15 for parietal cortex {y — —0.15a; -|- 0.01); (C) —0.15 for sensorimotor cortex 
{y = -0.15x + 0.02); (D) -0.15 for temporal cortex {y = -0.15x + 0.07). (E) Average 
glucose utilization rate of the entire cerebral cortex yields the specific exponent —0.15 
(y = -0.15X + 0.03). 

Fig. 3 

Scaling of the specific glucose utilization rate in subcortical gray matter with brain 
volume. The specific metabolic scaling exponent has the following values: (A) —0.15 
for thalamus {y = —0.15a; + 0.03); (B) —0.14 for hippocampus {y = —0.14a; — 0.13), 
which represents limbic structures; (C) —0.15 for caudate {y — — 0.15a; -|- 0.02), which 
represents basal gangha; (D) —0.15 for cerebellum (y = —0.15a; — 0.09). 
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Fig. 4 

Scaling of the volume-specific glucose utilization rate in white matter with brain volume. 

(A) Corpus callosum metabohsm yields the exponent —0.23 (y = —0.23x — 0.45), and 

(B) internal capsule has a similar exponent —0.24 (y = — 0.24x — 0.41). 
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